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Abstract 

Recent diffractive structure function measurements by the HI and ZEUS ex- 
periments at HERA are reviewed. Various data sets, obtained using system- 
1/-^ . atically different selection and reconstruction methods, are compared. NLO 

O I DGLAP QCD fits are performed to the most precise HI and ZEUS data and 

diffractive parton densities are obtained in each case. Differences between the 
dependences of the HI and ZEUS data are reflected as differences between 
Q i the diffractive gluon densities. 

1 Introduction 

In recent years, several new measurements of the semi-inclusive 'diffractive' deep inelastic scattering 
(DIS) cross section for the process ep eXY at HERA have been released by the HI and ZEUS 
^ ■ experiments [1-6]. The data are often presented in the form of a t-integrated reduced diffractive neutral 

^T) ■ current cross section aP^^\ defined through' 
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^ ■ or in terms of a diffractive structure function F^^^\x]p, f3, Q"^). Neglecting any contributions from 

^ . exchange, 

such that o"r = is a good approximation except at very large y. The new data span a wide 

kinematic range, covering several orders of magnitude inQ"^, (3 and xp. 

^ ' Within the framework of QCD hard scattering coUinear factorisation in diffractive DIS [7], these 

■ data provide important constraints on the diffractive parton distribution functions (dpdf 's) of the proton. 

These dpdf 's are a crucial input for calculations of the cross sections for less inclusive diffractive pro- 
cesses in DIS, such as dijet or charm production [8,9]. In contrast to the case of inclusive scattering, the 
dpdf 's extracted in DIS are not expected to be directly applicable to hadron-hadron scattering [7, 10-12]. 
Indeed, diffractive factorisation breaks down spectacularly when HERA dpdf 's are applied to diffrac- 
tive proton-proton interactions at the TEVATRON [13]. It may, however, be possible to recover good 
agreement by applying an additional 'rapidity gap survival probability' factor to account for secondary 
scattering between the beam remnants [14-17]. The HERA dpdf 's thus remain an essential ingredient in 
the prediction of diffractive cross sections at the LHC, notably the diffractive Higgs cross section [18]. 
Although the poorly known rapidity gap survival probability leads to the largest uncertainty in such cal- 
culations, the uncertainty due to the input dpdf 's also plays a significant role. In [3], the HI collaboration 
made a first attempt to assess the uncertainty from this source, propagating the experimental errors from 
the data points to the 'HI 2002 NLO fit' parton densities and assessing the theoretical uncertainties from 
various sources. 

In this contribution, we investigate the compatibility between various different measurements of 
by HI and ZEUS. We also apply the techniques developed in [3] to ZEUS data in order to explore 
the consequences of differences between the HI and ZEUS measurements in terms of dpdf 's. 

'Contributed to the Proceedings of tlie Worksliop on HERA and tlie LHC, DESY and CERN, 2004-2005. 
'For a full definition of all terms and variables used, see for example [3]. 



2 Diffractive Selection Methods and Data Sets Considered 

One of the biggest challenges in measuring diffractive cross sections, and often the source of large 
systematic uncertainties, is the separation of diffractive events in which the proton remains intact from 
non-diffractive events and from proton-dissociation processes in which the proton is excited to form a 
system with a large mass, M^. Three distinct methods have been employed by the HERA experiments, 
which select diffractive events of the type ep —>■ eXY, where y is a proton or at worst a low mass proton 
excitation. These methods are complimentary in that their systematics due to the rejection of proton 
dissociative and non diffractive contributions are almost independent of one another. They are explained 
in detail below. 

• Roman Pot Spectrometer Method. Protons scattered through very small angles are detected 
directly in detectors housed in 'Roman Pot' insertions to the beampipe well downstream the inter- 
action point. The proton 4-momentum at the interaction point is reconstructed from the position 
and slope of the tracks in these detectors, given a knowledge of the beam optics in the intervening 
region. The Roman Pot devices are known as the Leading Proton Spectrometer (LPS) in the case 
of ZEUS and the Forward Proton Spectrometer (EPS) in HI. The Roman pot method provides 
the cleanest separation between elastic, proton dissociative and non-diffractive events. However, 
acceptances are rather poor, such that statistical uncertainties are large in the data sets obtained so 
far. 

• Rapidity Gap Method. This method is used by HI for diffractive structure function measurements 
and by both HI and ZEUS for the investigation of final state observables. The outgoing proton is 
not observed, but the diffractive nature of the event is inferred from the presence of a large gap in 
the rapidity distribution of the final state hadrons, separating the X system from the unobserved Y 
system. The diffractive kinematics are reconstructed from the mass of the X system, which is well 
measured in the main detector components. The rapidity gap must span the acceptance regions of 
various forward^ detector components. Eor the HI data presented here, these detectors efficiently 
identify activity in the pseudorapidity range 3.3 < t] < 7.5. The presence of a gap extending 
to such large pseudorapidities is sufficient to ensure that ^1.6 GeV. In light of the poor 
knowledge of the spectrum at low masses, no attempt is made to correct the data for the small 
remaining proton dissociation contribution, but rather the cross sections are quoted integrated over 
My < 1.6 GeV. 

• Method. Again the outgoing proton is not observed, but rather than requiring a large rapidity 
gap, diffractive events ai^e selected on the basis of the inclusive In distribution. Diffractive 
events are responsible for a plateau in this distribution at low In , such that they can be se- 
lected cleanly for the lowest values. At intermediate M^, non-diffractive contributions are 
subtracted on the basis of a two component fit in which the non-diffractive component rises expo- 
nentially. This method is used for diffractive structure function measurements by ZEUS. It does 
not discriminate between elastic and low proton-dissociative contributions. Results are quoted 
for My < 2.3 GeV. 

Four recent data sets are considered, for which full details of luminosities and kinematic ranges can be 
found in tabled 

• Published data from ZEUS taken in 1998 and 1999, using the M^ method and taking advantage 
of the increased forward acceptance offered by a new plug calorimeter ('ZEUS-M^^') [1]. 

• Published ZEUS data obtained with the LPS using data taken in 1997 ('ZEUS-LPS') [2]. 

• Preliminary HI data obtained using the rapidity gap method, combining three measurements using 
different data sets from the period 1997-2000 for different regions in ('Hl-LRG') [3-5]. 

• Preliminary HI data obtained using the FPS, based on data taken in 1999 and 2000 ('Hl-FPS') [6]. 



^The forward hemisphere corresponds to that of the outgoing proton beam, where the pseudorapidity -q — — in tan 6*72 is 
positive. 
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ZEUS-A/x 


[1] 


My- method 


4.2 


< 2.3 


2.7... 55 


0.003... 0.975 


0.0001 . 


. .0.03 


ZEUS-LPS 


[2] 


Roman Pot 


12.8 


Mp 


2.4... 39 


0.007... 0.48 


0.0005 . 


. .0.06 


Hl-LRG 


[3-5] 


Rapidity Gap 


3.4... 63 


< 1.6 


1.5... 1600 


0.01... 0.9 


0.0001 . 


..0.05 


HI -EPS 


[6] 


Roman Pot 


25 


Mp 


2.6... 20 


0.01... 0.7 


0.002 . . 


.0.05 



Table 1: Overview of the data sets discussed here. The quoted kinematic ranges inQ^, f} and xp correspond to the bin centres. 
3 Comparisons between Data Sets 

In this section, the x jp dependences of the data from the different measurements are compared at fixed 
values of and /3. Since the various measurements are generally presented at different and /? 
values, it is necessary to transport the data to the same values. The (3 and values of the Hl-LRG 
data are chosen as the reference points. The factors applied to data points from the other measurements 
are evaluated using two different parameterisations, corresponding to the results of QCD fits to 1994 HI 
data [19] and to a subset of the present Hl-LRG data at intermediate [3] (see also sectionlU. In order 
to avoid any significant bias arising from this procedure, data points are only considered further here if 
the correction applied is smaller than 50% in total and if the correction factors obtained from the two 
parameterisations are in agreement to better than 25%. In practice, these criteria only lead to the rejection 
of data points in the ZEUS-Mj. data set at = 55 GeV^ and /? = 0.975, where the poorly known high /? 
dependence of the diffractive cross section implies a large uncertainty on the factors required to transport 
them to (3 = 0.9. Elsewhere, there is reasonable agreement between the factors obtained from the two 
parameterisations and no additional uncertainties are assigned as a consequence of this procedure. 

Since the various data sets correspond to different ranges in the outgoing proton system mass, 
My , additional factors are required before comparisons can be made. For all data and fit comparisons, 
all data are transported to the HI measurement range of My < 1.6 GeV and |t| < 1 GeV^. The leading 
proton data are scaled by a factor 1.1 [20] to correspond to the range < 1.6 GeV and the ZEUS- 
Mj, data are scaled to the same range by a further factor of 0.7 [1], such that the overall factor is 0.77. 
The uncertainties on these factors are large, giving rise to normalisation uncertainties of perhaps 15% 
between the different data sets. 

The ZEUS-LPS and Hl-FPS data are compared in figure^ Within the experimental uncertainties, 
the two data sets are in good agreement. Both data sets are also consistent with a parameterisation of 
the Hl-LRG data [3] based on the HI 2002 NLO QCD fit, which is also shown. This good agreement 
between the Hl-LRG and the Roman Pot data is also shown explicitly in figure |3l 

In figure|2l a comparison is made between the Hl-LRG and the ZEUS-Af,^ data after all factors 
have been applied. For much of the kinematic range, there is tolerable agreement between the two data 
sets. However, there are clear regions of disagreement. One is at the largest (3 (smallest M^), where the 
HI data lie significantly above the ZEUS data for < 20 GeV^. Another is at intermediate and low f3, 
where the two data sets show significantly different dependences on Q^. With the factor of 0.77 applied 
to the ZEUS data, there is good agreement at low Q^, but the ZEUS data lie below the HI data at large 
Q^. If the factor of 0.77 is replaced with a value closer to unity, the agreement improves at lai^ge Q^, but 
the HI data lie above the ZEUS data at low Q^. These inconsistencies between the different data sets are 
discussed further in section ID 

For completeness, figure |3l shows a comparison between all four data sets considered. 
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Fig. 1: Comparison of tlie Roman Pot data from HI and ZEUS, scaled by a factor 1.1 such that they correspond to < 
1.6 GeV. The and (3 values have been shifted to the Hl-LRG bin centres using small translation factors. The upper and 
lower curves form an error band on the predictions from the HI 2002 NLO QCD fit to the Hl-LRG data (experimental errors 
only). Dotted lines are used for kinematic regions which were not included in the fit. Normalisation uncertainties of ^Jq% on 
the ZEUS LPS data and 15% on the factor applied to shift the datasets to My < 1.6 GeV are not shown. 



4 Diffractive Parton Distributions 

4.1 Theoretical Framework and Fit to Hl-LRG Data 

In this contribution, we adopt the fitting procedure used by HI in [3], where next-to-leading order (NLO) 
QCD fits are performed to diffractive reduced cross section, a!^^^\ data [3,21] with 6.5 < < 
800 GeV^ and the /? and xjp ranges given in tabled 

The proof that QCD hard scattering coUinear factorisation can be applied to diffractive DIS [7] 
implies that in the leading log(Q^) approximation, the cross section for the diffractive process ep — 
eXY can be written in terms of convolutions of universal partonic cross sections cj^* with diffractive 
parton distribution functions (dpdf's) fP [11,22,23], representing probability distributions for a parton 
i in the proton under the constraint that the proton is scattered with a particular 4 momentum. Thus, at 
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Fig. 2: Comparison of tlie ZEUS-Mx data with a subset of the Hl-LRG data. The and /3 values of the ZEUS data have been 
shifted to the HI bin centres using small translation factors. The ZEUS data have also been multiplied by a universal factor of 
0.77, such that both data sets correspond to Afy < 1.6 GeV. Normalisation uncertainties of 15% on this factor and of ±6.7% 
on the HI data are not shown. 



leading twist,^ 



d^a{x,Q'^,xip,tYP^^^P' 
dxjp dt 



E 



(3) 



This factorisation formula is valid for sufficiently large and fixed Xip and t. It also applies to the case 
of proton dissociation into a system of fixed mass and thus to any cross section which is integrated 
over a fixed range in My . The partonic cross sections a'^^ are the same as those for inclusive DIS and the 
dpdf 's fP, which are not known from first principles, should obey the DGLAP evolution equations [25]. 

In addition to the rigorous theoretical prescription represented by equation an additional as- 
sumption is necessary for the HI fits in [3], that the shape of the dpdf's is independent of xjp and t 
and that their normalisation is controlled by Regge asymptotics [26]. Although this assumption has no 

''a framework also exists to include higher order operators [24]. 
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Fig. 3: Summary plot of all diffractive DIS data sets considered here. Additional Hl-LRG data with < 2.5 GeV^, 

2 2 2 2 2 

Q =45 GeV and Q > 60 GeV are not shown. The Q and /? values for all data sets have been shifted to the HI bin 
centres using small translation factors. The ZEUS data have been multiplied by a universal factor of 0.77 and the LPS and FPS 
data by factors of 1.1, such that all data sets correspond to My < 1.6 GeV. Relative normalisation uncertainties of 15% due 
to these factors and further normalisation uncertainties of ±6.7% (Hl-LRG) and II^qS (ZEUS-LPS) data are not shown. 



solid basis in QCD, it is compatible with the data fitted. The diffractive parton distributions can then be 
factorised into a term depending only on xjp and t and a term depending only on x (or (3) and Q^: 



ffixjp, t, X, Q^) = fjp/p{xjp, t) ■ ff{(3 = x/xip, g2) . (4) 

Under this 'Regge' factorisation assumption, the diffractive exchange can be treated as an object (a 
'pomeron', IP) with a partonic structure given by parton distributions ff (/3, Q^). The variable /? then 
corresponds to the fraction of the pomeron longitudinal momentum carried by the struck parton. The 
'pomeron flux factor' f]p/p{x]p,t) represents the probability that a pomeron with particular values of 
xjp and t couples to the proton. 

In the fit, the x jp dependence is parameterised using a Regge flux factor 

ftmin e^^* 

fs^/pix]p, t)=A- j 2ar(t)-l ' 

cut X JP 

where tout = —10 GeV^, |imm| is the minimum kinematically allowed value of |t| and the pomeron 
trajectory is assumed to be linear, aip{t) = aip{0) + a'jpt. The parameters Bp and a' and their uncer- 
tainties are fixed as described in [3]. The value of A is chosen such that the flux factor is normalised to 
unity atx]p = 0.003. The pomeron intercept is then obtained from the xjp dependence of the data and 
takes the value a^{0) = 1.173 ± 0.018 (stat.) ± 0.017 (syst.) (model). 

The description of the data is improved with the inclusion of an additional separately factorisable 

sub-leading exchange with a trajectory intercept of ajR(O) = 0.50 and parton densities taken from a 
parameterisation of the pion [27]. This exchange contributes significantly only at low P and large xjp. 

The dpdf 's are modelled in terms of a light flavour singlet 

E(z) = u{z) + d{z) + s{z) + u{z) + d{z) + s{z) , (6) 

with u = d = s = u = d = s and a gluon distribution g{z) at a starting scale Qq = 3 GeV^. Here, 
z is the momentum fraction of the parton entering the hard sub-process with respect to the diffractive 
exchange, such that z = P for the lowest-order quark parton model process, whereas < /3 < 2; for 
higher order processes. The singlet quark and gluon distributions are parameterised using the form 



Y,C}Pj{2z-l) 



2 



0.01 

e—^ , (7) 



where Pj{0 is the j^^ member of a set of Chebychev polynomials^. The series is squared to ensure 
positivity. The exponential term is added to guarantee that the dpdf 's tend to zero in the limit of z ^ 1. 
It has negligible influence on the extracted partons at low to moderate z. The numbers of terms in the 
polynomial parameterisations are optimised to the precision of the data, with the first three terms in the 
series used for both the quark singlet and the gluon distributions, yielding 3 free parameters (Cf and 
Cj) for each. The normalisation of the sub-leading exchange contribution at high xjp is also determined 
by the fit such that the total number of free parameters is 7. The data used in the fit are restricted to 
> 2 GeV to suppress non-leading twist contributions. The effects of F[' are considered through its 
relation to the NLO gluon density, such that no explicit cut on y is required. 

The NLO DGLAP equations are used to evolve the dpdfs to > Qq using the method of [28], 

extended for diffraction. No momentum sum rule is imposed. Charm quarks are treated in the massive 
scheme (appearing via boson gluon fusion processes) with rric = 1.5 ± 0.1 GcV. The strong coupling is 
set via^ ^qcd — 200 ± 30 MeV. The statistical and experimental systematic errors on the data points 

^Pi = 1, P2 = 5 and P,+i(0 = 2^P,{0 - Pi-i(C). 

^Although this value is rather different from the world average, we retain it here for consistency with previous HI prelimi- 
nary results, where it has been used consistently for QCD fits [3] and final state comparisons [8, 9]. 



and their correlations are propagated to obtain error bands for the resulting dpdfs, which correspond to 
increases in the by one unit [29]. A theoretical error on the dpdfs is estimated by variations of Aqcd> 
nic and the parameterisation of the xp dependences as described in [3]. No theoretical uncertainty is 
assigned for the choice of parton parameterisation, though the results are consistent within the quoted 
uncertainties if alternative approaches [30] are used. No inhomogeneous term of the type included in [31] 
is considered here. The presence of such a term would lead to a reduction in the gluon density extracted. 

The central fit gives a good description of the data, with a of 308.7 for 306 degrees of freedom. 
The resulting diffractive quark singlet and gluon distributions are shown in figure^] Both extend to large 
fractional momenta z. Whereas the singlet distribution is well constrained by the fit, there is a substantial 
uncertainty in the gluon distribution, particularly for z > 0.5. The fraction of the exchanged momentum 
carried by gluons integrated over the range 0.01 < z < 1 is 75 ± 15% (total error), confirming the 
conclusion from earlier work [19] that diffraction is a gluon-induced phenomenon. These dpdfs have 
been astonishingly successful in describing diffractive final state data in DIS such as charm [9] and jet [8] 
production, which, being induced by boson-gluon fusion-type processes, are roughly proportional to the 
diffractive gluon density. 

4.2 Fit to ZEUS Data 

A very similar fit to that described in section l4!n is performed to the ZEUS-Mj, data and the implications 
of the differences between the data sets to the dpdfs are investigated. The data are fitted in their original 
binning scheme, but are scaled to My < 1.6 GeV using the factor of 0.77. As for the fit to the HI 
data, the first 3 terms are included in the polynomial expansions for the quark and gluon densities at the 
starting scale for QCD evolution. The same fit program, prescription and parameters are used as was the 
case for the HI 2002 NLO fit, with the following exceptions. 

• ZEUS-Mx data with > 4 GeV^ are included in the fit, whereas only HI data with > 
6.5 GeV^ are included. It has been checked that the result for ZEUS is not altered significantly if 
the minimum value is increased to 6 GeV^. 

• The quadratic sum of the statistical and systematic error is considered, i.e. there is no treatment of 
correlations between the data points through the systematics. 

• No sub-leading Reggeon exchange component is included in the parameterisation. Including one 
does not improve or alter the fit significantly. 

• The Pomeron intercept is fitted together with the dpdfs, in contrast to the two stage process of [3]. 
This does not influence the results significantly, though it does decrease the uncertainty on a-^ (0). 

The fit describes the ZEUS-M^ data well (x^ = 90 for 131 de grees of freedom) and yields a value for 
the Pomeron intercept of ajp(O) = 1.132 it 0.006 (experimental error only). This value is in agreement 
with the HI result if the full experimental and theoretical errors are taken into account. A good fit is thus 
obtained without any variation of ajp(O) with Q"^ or other deviation from Regge factorisation. 

The diffractive parton densities from the fit to the ZEUS-M^^^ data are compared with the results 
from HI in figure |4l The differences observed between the HI and the ZEUS data are directly reflected 
in the parton densities. The quark singlet densities are closely related to the measurements of them- 
selves. They are similar at low where the HI and ZEUS data are in good agreement, but become 
different at larger Q^, where discrepancies between the two data sets are observed. This difference be- 
tween the dependences of the HI and ZEUS data is further reflected in a difference of around a 
factor of 2 between the gluon densities, which are roughly proportional to the logarithmic derivative 
dFi^/dlnQ'^ [32]. 

The Hl-LRG and ZEUS-M^ data are shown together with the results from both QCD fits in figure 
|5] Both fits give good descriptions of the data from which they are obtained. The differences between 
the two data sets are clearly reflected in the fit predictions, most notably in the dependence. 



NLO QCD fits to H1 and ZEUS data 




z z 

H1 2002 NLO fit (prel.) 
(exp. error) 
(exp.+theor. error) 

NLO fit to ZEUS-IVIx (exp. error) 

Fig. 4: Diffractive quark singlet and gluon pdf's for various values, as obtained from the NLO DGLAP fits to the Hl-LRG 
and ZEUS-iWx data. The bands around the HI result indicate the experimental and theoretical uncertainties. The dotted lines 
around the result for ZEUS indicate the experimental uncertainty. The ZEUS data used in the fit are scaled by a normalisation 
factor of 0.77 to match the Hl-LRG range of < 1.6 GeV. This factor is reflected in the normalisations of the quark and 
gluon densities. An uncertainty of 15% on this factor is not included in the error bands shown. 

5 Summary 

Recent diffractive structure function data from HI and ZEUS have been compared directly. The leading 
proton data from both experiments (Hl-FPS and ZEUS -LPS) are in good agreement with one other and 
with the HI large rapidity gap data (Hl-LRG). There is reasonable agreement between the Hl-LRG and 
the ZEUS-Mj, data over much of the kinematic range. However, differences are observed at the highest 
P (smallest M^) and the dependence at intermediate to low /5 is weaker for the ZEUS-M^, data than 
is the case for the Hl-LRG data. 

An NLO DGLAP QCD fit has been performed to the ZEUS-M,^ data, using the same theoreti- 
cal framework, assumptions and parameterisations as have been employed previously for the HI -2002- 
prelim NLO QCD fit to a subset of the Hl-LRG data. As a consequence of the differences between the 
dependences of the Hl-LRG and ZEUS-M-^ data, the gluon density obtained from the ZEUS data is 
significantly smaller than that for HI. 
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